First-principles study of the electrooptic effect in ferroelectric oxides 
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Wc first present a method to compute the electrooptic tensor from first principles, explicitly tak- 
ing into account the electronic, ionic and piezoelectric contributions. We then study the non-linear 
optic behavior of three paradigmatic ferroelectric oxides. Our calculations reveal the dominant con- 
tribution of the soft mode to the electrooptic coefficients in LiNbOs and BaTiOa and its minor role 
in PbTiOa. We identify the coupling between the electric field and the polar atomic displacements 
along the B-0 chains as the origin of the large electrooptic response in perovskite ABO3 compounds. 
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The electrooptic (EO) effect describes the change of 
refractive index of a material in a static electric field and 
is exploited in various technological applications pj . Fer- 
roelectric ABO3 compounds exhibit unusually large EO 
coefficients and are therefore materials of choice for op- 
tical devices. Since the seventies, LiNb03 EO modula- 
tors have been widely used in fiber-optic transmission 
systems 0- More recently, there has been increasing in- 
terest in epitaxially grown BaTiOa thin films for optical 
waveguide modulators The EO effect is the origin of 
the photorefractive effect, exploited in non- volatile holo- 
graphic data storage in LiNbOs Q. It was also used to 
probe locally the ferroelectricity in SrTi03 films 

Finding better EO materials is a desirable goal. How- 
ever, the experimental characterization of optical non- 
linearities requires high-quality single crystals that are 
not always directly accessible nor easy to make. Input 
from accurate theoretical calculations allowing to predict 
the non-linear optical behavior of crystalline solids would 
therefore be particularly useful. 

For many years, theoretical investigations of non-linear 
optical phenomena were restricted to semi-empirical ap- 
proaches such as shell models or bond-charge mod- 
els In the last decade, significant theoretical 
advances have been reported concerning first-principles 
density functional theory (DFT) calculations of the be- 
havior of periodic systems in an external electric field 0, 
ITo| and opened the way to direct predictions of various 
optical phenomena. Recently, particular attention has 
been paid to the calculation of non-linear optical (NLO) 
susceptibilities and Raman cross sections 0, ^| . 

In this Letter, we go one step further and present a 
method to predict the linear EO coefficients of periodic 
solids within DFT. Our method is very general, and can 
be applied to paradigmatic ferroelectric oxides : LiNbOs, 
BaTi03 and PbTi03. We find that first-principles calcu- 
lations are fully predictive, and provide significant new 
insights into the microscopic origin of the EO effect. In 
particular, we highlight the predominent role of the soft 
mode in the EO coupling of LiNb03 and BaTiOs, in con- 
trast with its minor role in PbTi03. 

At linear order, the dependence of the optical dielectric 
tensor e.^ on the static (or low-frequency) electric field 



£ 1 is described by the linear EO tensor r 
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Throughout this paper, we follow the convention of us- 
ing Greek and Roman indexes (resp.) to label static and 
optical fields (resp.). We write all vector and tensor com- 
ponents in the system of cartesian coodinates defined by 
the principal axes of the crystal under zero field. We also 
refer to the atomic displacements r KQ [k labels an atom 
and a a cartesian direction] within the basis defined by 
the zone-center transverse optic (TO) phonon eigendis- 
placements u m (na): r Ka = J2,n T mU m (na). 

Let us first consider the clamped (zero strain) EO ten- 
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sor, r , in which all electric-field induced macroscopic 



strains 77 are forbidden. This is achieved experimentally 
by working at a frequency sufficiently high to avoid strain 
relaxations but low compared to the frequency of the TO 
modes. Within the Born-Oppenheimer approximation, 
we express the total derivative of as the sum of two 
partial derivatives with respect to E 1 and r m : 



d£~, 



d£^ 



T = 



dy {1) 

u 1 rn 



£=0 



d£~, 



(2) 



The derivative in the first term of the righthand side of 
Eq. is computed at clamped atomic positions. It de- 
scribes the electronic contribution to the EO tensor and 
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is proportional to the NLO susceptibilities Xiji 
ond term represents the ionic contribution. It depends 
on the first-order change of the linear dielectric suscep- 
tibility due to atomic displacements, and is related to 

the Raman susceptibility a™ = V^J2 KOl 7fF 
mode m [Q is the unit cell volume], as well as to 
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amplitude of the ionic relaxation induced by the field 
£y. dr m /d£ 7 can be expressed in terms of (i) the TO 
phonon mode frequencies uj m and (ii) the TO mode po- 
larities Pm.-y — Y^k ^k, <yfi^m( K /3)> directly linked to the 
infrared (IR) intensities [13. Combining this with the 
previous equations, we obtain the clamped EO tensor 0] 
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where and rij are the principal refractive indices. 

Let us now consider the undamped (zero stress) EO 
tensor, rL . It can be shown 1141 that the macroscopic 
expression proposed in Ref. [l5j is still valid at the mi- 
croscopic level so that the the piezoelectric contribution 
to rfj can be computed from the elasto-optic coefficients 
Vijap and the piezoelectric strain coefficients <l ia p 
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An expression similar to Eq. (j3J was previously used 
by Johnston 0] to estimate the clamped EO tensor of 
LiNb03 and LiTa03 from IR and Raman measurements. 
However, this semi-empirical approach was limited by the 
indeterminacy of the relative sign of p m and ct m . As 
discussed below, the direct evaluation of Eq. (0 and (@J 
from first principles provides an easier and more accurate 
estimate of the EO tensor. 

We have implemented this formalism in the ABINIT 
open software |17j . within the local density approxima- 
tion (LDA) to the DFT. The optical dielectric tensor, 
Born effective charges, phonon frequencies and eigendis- 
placements are computed from linear response |13| . The 
piezoelectric strain coefficients d 7Ct ^ are deduced from 
the piezoelectric stress coefficients e 7Ct ^ and the elastic 
constants. These two quantities, as well as the elasto- 
optic tensor Pij a /3, are obtained from finite differences. 

The non-linear response functions xtji an d / dr Ka 
are computed from a perturbative approach using a new 
implementation based on the 2n + 1 theorem. To reach 
reasonable k-point sampling convergence, we combined 
the recently proposed PEAD expression [9j and the fi- 
nite difference formula of Marzari and Vanderbilt |l8| to 
compute the perturbation expansion of the polarization. 
More details will be provided elsewhere 0] . The method 
was tested on various cubic semiconductors and provides 
results in close agreement with earlier studies [llj, uM ■ 

For BaTiOs and PbTiOs, we use extended norm- 
conserving pseudopotentials |l9j | , a planewave kinetic en- 
ergy cutoff of 45 hartree and a 10 x 10 x 10 k-point grid. 
For LiNbOa , we use the same norm-conserving pseudopo- 
tentials as in Ref. [2£j as well as the Born effective 
charges, phonon frequencies and eigenvectors already re- 
ported in that paper. For this compound, a 8 x 8 x 8 
k-point grid and a planewave kinetic energy cutoff of 35 
hartree give converged values for an d dytij ldr KOL . 

First, we study LiNb03. This compound has a trigonal 
symmetry with 10 atoms per unit cell. The theoretical 
lattice constants and atomic positions are reported in 
Ref. [2(j. LiNb03 undergoes a single transition at 1480 
K from a centrosymmetric high-T paraelectric Ric phase 



to a ferroelectric low-T R3c ground state. The form of the 
EO tensor depends on the choice of the cartesian axes. 
Here, we follow the I.R.E. Piezoelectric Standards plj . 

With this choice of axes, the EO tensor in the ferroelec- 
tric phase of LiNb03 has 4 independent elements (Voigt 
notations): ri3, r33, r 22 and rsi. The TO modes can be 
classified into AA\ + 5A 2 + 9E. The Ai and E modes are 
Raman and IR active. Only the Ai modes couple to ri3 
and r33, while the E modes are linked to r 22 and rsi. Ta- 
ble [I] gives these four clamped coefficients |22 . as well as 
the contribution of each optical phonon. For comparison, 
we mention the coefficients computed by Johnston [l6| 
from measurements of IR and Raman intensities (IR + 
R) as well as the results of a bond-charge model (BCM) 
calculation by Shih and Yariv The first-principles 
calculations correctly predict the sign of the four EO co- 
efficients [2l]]. The absolute values are also well repro- 
duced by our method, especially if we take into account 
the fact that NLO properties are generally difficult to 
determine accurately. The experimental values are sensi- 
tive to external parameters such as temperature changes 
[23| and the stoichiometry of the samples. For exam- 
ple, using crystals of various compositions, Abdi and co- 
workers measured absolute values between 1.5 pm/V and 
9.9 pm/V for r 22 0- These difficulties support the need 
for sophisticated theoretical tools to predict NLO prop- 
erties. In contrast to the models of Refs. || our 
method is predictive and does not use any experimen- 
tal parameters. Moreover, it reproduces r^ 3 , r 33 and r 22 
better than the semiempirical models. 

Our approach also provides some insight into the ori- 
gin of the high LiNb03 EO response. All EO coefficients 
are dominated by the ionic contribution of the Ai TOl 
and the E TOl modes. This can be explained as fol- 
lows. At the paraelectric-ferroelectric phase transition, 
the unstable A 2u and E„ modes of the paraelectric phase 
transform to low-frequency and highly polar modes in the 
ferroelectric phase |2£j, generating a large EO response 
if they exhibit, in addition, a large Raman susceptibil- 
ity. The Ai TOl and E TOl modes of the ferroelectric 
phase have a strong overlap of respectively 0.82 and 0.68 
with the unstable A 2u and E u modes of the paraelectric 
phase and combine giant polarity |2(i| and large Raman 
susceptibility (see below for the Ai mode). 

In Table |U we also report the undamped EO coeffi- 
cients in LiNb03. As the piezoelectric coefficients d3i 
(-1 pC/N) and d33 (6 pC/N) are small compared to dis 
(55.9 pC/N) and d 22 (21.6 pC/N), the piezoelectric ef- 
fect is important for r 22 and rj^ and negligible for rf 3 and 
rfg. The undamped EO coefficient is nearly twice 
as large as the clamped one. Moreover, its theoretical 
value is in better agreement with the experiment than 
that of the clamped one. This suggests that the piezo- 
electric contribution was not entirely eliminated during 
the measurement of r^ ; the correct value of the clamped 
coefficient might be closer to the theoretical 14.9 pm/V. 

Second, we study PbTi03 and BaTiOs. Both com- 
pounds are stable at room temperature in a ferroelec- 
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TABLE I: EO tensor (pm/V) in LiNb0 3 : electronic, ionic 
and piezoelectric contributions, and comparison with experi- 
ment, for the clamped and undamped cases. The ionic part 
is split into contributions from TO modes (u> m in cm -1 ). 
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trie distorted perovskite structure of tetragonal P4mm 
symmetry with 5 atoms per unit cell |26| . In the P4mm 
phase, the TO modes can be classified into 3Ai+4E+Bi. 
The EO tensor has only three independent elements: ri3, 
and r33, coupling to the Ai modes, and r42, linked to the 
E modes. The Bi-mode is IR inactive and does not in- 
fluence the EO tensor. The results are shown in Table ITT1 
For PbTiOa, we found only measurements of r^ 3 and 
r 33 , which agree well with our theoretical results. More- 
over, our calculation predicts that PbTi03 exhibits a 
large tJ 2 , in spite of its low r 33 . Combined with other 
advantageous features, such as small thermo-optic coeffi- 
cients [29j . this suggests that PbTiOs might be an inter- 
esting candidate for EO applications if properly oriented. 

In BaTiOs, the low temperature structure is rhombo- 
hedral. The P4mm phase is unstable and exhibits, in the 
harmonic approximation, an unstable E-mode that pre- 
vents the use of Eq. to compute x\ 2 . The theoretical 
estimates of r^ 3 and r 33 are reasonably accurate despite 
an underestimation of the theoretical r 33 . The origin of 
the error can be attributed to various sources. First, 
the values computed for the P4mm phase correspond to 
an extrapolation of the EO tensor to K, while experi- 
mental results are obtained at room temperature. Also, 
linear and NLO susceptibilities can be relatively inaccu- 
rate within the LDA. In this context, note the use of the 
LDA optical refractive indexes in Eq. J3J), overestimating 
the experimental values by about 10 %. 

We compare now the NLO response of the three com- 
pounds. r^ 3 is similar for all of them, while r 33 is signif- 
icantly smaller in PbTiOs than in LiNb03 and BaTiOs. 



TABLE II: Electronic and ionic contributions of individual 
TO modes (uj m in cm" 1 ) to the clamped EO tensor (pm/V) 
in the P4mm phase of PbTi0 3 and BaTi0 3 . 

PbTiOs BaTiOs 
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TABLE III: Raman susceptibilities and mode polarities (10 2 
a. u.) of the Ai TO modes in LiNb0 3 , BaTi0 3 and PbTi0 3 . 

LiNbOs BaTiOs PbTiOs 

P3 Ctll Q33 ps Qll Q33 p-j an Q33 

TOl 3.65 -0.70 -2.02 1.22 -0.16 -0.13 1.25 -0.67 -0.43 

T02 0.45 0.30 0.53 3.25 -1.18 -2.73 2.18 -0.75 -0.33 

T03 0.67 0.18 -0.05 1.74 -1.26 -2.55 2.68 -2.42 -2.28 

T04 3.82 -1.96 -3.23 



In the latter two compounds, the magnitude of r 33 is 
dominated by one particular phonon mode. In BaTiOs, 
the T02 mode at 300 cm -1 has a similar strong overlap 
(92%) with the unstable mode in the paraelectric phase 
than the TOl modes in LiNbOs, as previously discussed. 
In PbTiOs, all Ai modes contribute almost equally to 
r 33 . The T02 mode at 357 cm" 1 has the strongest over- 
lap (73%) with the soft mode in the cubic phase. Surpris- 
ingly, its contribution to r 33 is 23.5 times smaller than 
the contribution of the T02 mode in BaTiOs. 

To identify the origin of the distinctive behavior of 
PbTiOs, we report in Table IIIII the mode polarities 
and Raman susceptibilities of the Ai TO modes. In 
the three compounds, a has two independent elements 
a ii and 033 that determine the amplitude of r^ 3 and 
r 33 . «33 is large for the TOl mode in LiNbOs and 
the T02 mode in BaTiOs. On the other hand, it is 
the smallest for the T02 mode in PbTi03, in agree- 
ment with experiments |30j. Combined with a higher 
frequency (^% bTi o 3 /^BaTi0 3 = 1A1 )> a lower polarity 
(PBaTi0 3 /ppbTi0 3 = 1-49), and a larger value of the re- 
fractive index (n>pbTi0 3 / n BaTi0 3 ~ 1-35), this weak Ra- 
man susceptibility (ctBaTi0 3 l otpvno 3 — 8.27) explains 
the weak contribution of the T02 mode to r 33 in PbTiOs. 

The microscopic origin of the lower Ai T02 mode Ra- 
man susceptibility in PbTiOs, compared to BaTiOs, is 
explained by the decomposition of 1233 into contributions 
of the individual atoms in the unit cell (see Table llVjl . 
In both perovskites, the major contributions to the Ra- 
man susceptibility of the Ai T02 modes are 033 (Ti) and 



4 



TABLE IV: Decomposition of the Raman susceptibility of 
the Ai T02 mode in BaTiOs and PbTiOa into contributions 
from the individual atoms in the unit cell [26| . 
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Q!33(Oi) |2^; ^33 is mostly due to the atomic displace- 
ments of the atoms located on the Ti-0 chains oriented 
along the polar direction. First, the derivatives of X33 
versus atomic displacement are of opposite sign for Ti 
and Oi atoms, and significantly larger in BaTiOa than 
in PbTi03. Second, the opposing displacements of Ti 
and Oi atoms in the T02 mode in BaTiOs produce con- 
tributions that add to yield a giant 033. On the other 
hand, the in-phase displacements of Ti and Oi in PbTiC>3 
produce contributions that cancel out, giving a small 033. 
This distinct behaviour goes beyond a simple mass effect. 



Changing the mass of Pb to that of Ba in the dynamical 
matrix of PbTiC>3 has no significant effect on the relative 
Ti— O displacement. Large atomic displacements of oppo- 
site direction along the Ti-0 chains are therefore needed 
to generate a large 033 and potentially a large r^s. 

In summary, we presented a method to compute the 
EO tensor from first principles. In LiNbC>3 and BaTiC>3, 
the large EO response originates in the giant contribu- 
tion of the successor of the soft mode, which combines low 
frequency, high polarity and high Raman susceptibility. 
In comparison, the contribution of the similar mode in 
tetragonal PbTiOs is rather weak due to its low Raman 
susceptibility. In the perovskites, the Raman susceptibil- 
ity is principally determined by the atomic displacements 
along the B-0 chains in the polar direction. This sug- 
gests that the search for new perovskite oxides with good 
EO properties should focus on compounds with large rel- 
ative B-0 atomic displacements along the chains. 
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